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cr> ■ Abstract 

> : 

The B-meson light-cone wavefunctions, ^>±(uj, z 2 ), are investigated up to the next-to-leading 

l> ■ 

order in Fock state expansion in the heavy quark limit. In order to know the transverse momentum 

t> : 

dependence of the B-meson wavefunctions with 3-particle Fock states' contributions, we make 
use of the relations between 2- and 3- particle wavefunctions derived from the QCD equations of 

^ : 

Qh! motion and the heavy quark symmetry, especially two constraints derived from the gauge field 

equation of motion are employed. Our results show that the use of gluon equation of motion can 
give a constraint on the transverse momentum dependence x(u>, of the B-meson wavefunctions, 

^ ■ 

whose distribution tends to be a hyperbola- like curve under the condition < c\ < 1, which 
is quite different from the WW type wavefunctions, whose transverse momentum dependence 
X^(w,kj_) is merely a delta function. Based on the derived results, we propose a simple model 
for the B-meson wavefunctions with 3-particle Fock states' contributions. 
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I. INTRODUCTION 



The non-perturbative light-cone (LC) wavefunction (WF) /distribtuion amplitude (DA) 
of the B meson plays an important role for making reliable predictions for exclusive B meson 



iMflflflflfla. 



decays. The B-meson DA has been investigated in various approaches 

Recently, Ref.fiol] claims that it is the B-meson WF rather than the B-meson DA that 
is more relevant to the B decays, and in the framework of the fc^-factorization theorem 



they proved that the B-meson WF is renormalizable after taking into account the 
renormalizat ion-group (RG) evolution effects, meanwhile, the undesirable feature JjJ of the 
B-meson DA under evolution can be removed. Since the B-meson WF still poses a major 
source of uncertainty in the study of the B decays, hence, theoretically, it is an important 
issue to study on it. 

Ref. Q, as well [kJ, 

presents an analytic solution for the B-meson WFs ty±(u, z 2 ), 
which satisfies the constraints from the QCD equations of motion and the h eavy - quark 

Qrf 
. It is found that in the Wandzura-Wilczek(WW) approximation [151] . which 

corresponds to the valence quark distribution, the B-meson WFs can be determined uniquely 

in terms of the "effective mass" , the A, defined in the Heavy Quark Effective Theory (HQET) 

. Ref.[17] shows that when taking A G (O.BGeV, 0.6GeV), one can give a reasonable 
perturbative QCD result for B —>■ 71 transition form factor that is consistent with what was 
obtained in the LC sum rule calculation Q| and the lattice QCD simulation [rgj ]. 

It should be noted that in the WW approximation, the obtained analytic results for 
the B-meson WFs ^±(uj,z 2 ) are unique, and the only missing part for practical numerical 
use is the RG evolution effect. However, there is very limited knowledge on the higher Fock 
states' contributions. In Ref.Q, the B-meson distributions <fi±(uj) with 3-particle Fock states 
are given and a rough estimation presented there shows that the 3-particle Fock states' 
contributions might considerably broaden the transverse momentum distribution that is 
derived from the WW approximation. Recently, based on the QCD sum rule analysis and 
taking only the two 3-particle distributions ^a(p, and \EV(P> consideration, Ref. [3] 
connects the asymptotic behavior of their difference to the well-known quark-antiquark- 
gluon DA (p3 n (a.i) {oti(i = 1,2,3) are the fractions of the pion momentum carried by the 
corresponding partons and satisfy «i + a<i + — 1) ; i-e. in the small p and £ region, 
*v(p,0 - Va(p,€) ~ V?37r(ai, 1 - ax - a3,a 3 )L = p/ mSjQ 3 = £/ ms ~ p£ 2 - 
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In this paper, we are going to investigate the B-meson WFs ^/±(u,z 2 ) with the contri- 
butions from 3-particle Fock states. From the heavy-quark symmetry and the equations of 
motion for the light degrees of freedom, we get several constraints on the behaviors of WFs. 
By adopting some assumptions, we will solve the B-meson WFs approximately, especially, 
by taking the aforementioned asymptotic behavior of the difference between ^a(p, and 
\EV(P) m t° account, we try to derive an explicit form for the B-meson distributions (j)±{oS) 
that include the 3-particle Fock states' contributions. 

The paper is organized as follows. In Sec. II, several differential equations for the B-meson 
WFs are obtained by using the heavy quark symmetry and the equations of motion for the 
light degrees of freedom. In Sec. Ill, approximate solutions for the B-meson WFs including 
the contributions from the 3-particle Fock states are investigated under three assumptions. A 
new model for the B-meson WFs and some discussions over its phenomenological implication 
are presented in Sec. IV. The last section is reserved for a summary. 



II. DIFFERENTIAL EQUATIONS FOR THE B-MESON WAVEFUNCTIONS 

In HQET, the B-meson WFs $?+(t, z 2 ) can be defined in terms of the vacuum-to-meson 

n 

matrix element of the nonlocal operators [20J: 

(i) 



(0\q( z )rh v (0)\B(p)) = --^Tr 



75 rHk(t,,V^ (MV ^ (< ' z2 



2 [ TW ' r 2t 
where z^ = (0, z~ , z±), z 2 = — z^_, v 2 = 1, t = v ■ z, and = Mt> M is the 4-momentum of 
the B meson with mass M, h v (x) denotes the effective 6-quark field and T is a generic Dirac 
matrix. The path-ordered gauge factors are implied between the constituent fields. Note 
that in the above definition, the separation between quark and antiquark is not restricted 
on the LC (z 2 = 0). For a fast moving meson, t — > oo, Eq. (JTJ) shows that ^/ + (t,z 2 ) is 
the leading-twist WF and z 2 ) is the subleading one. To know more about the twist 

structures for the B meson, readers are recommended to refer to Ref.(f| for details, where 
the relation between the geometric twist and the dynamic twist was discussed. 

In the heavy quark limit, the general Lorentz decomposition of the 3-particle matrix 
elements can attribute to four independent 3-particle WFs similar to the 3-particle LC DAs 



12j , i.e. i& v (t,u, z 2 ), ^ A (t,u,z 2 ), X A (t,u, z 2 ) and Y A (t, u, z 2 ): 
(O|?0*) gG^iuz) zTh v (0)\B(p)) 



1 



75 r 



1+ 



(2) 



■ a^z" ty v (t,u,z 2 



z\ 
t 



X A (t,u,z 2 ) +[tt f 



t 



Y A (t,u,z 2 ) 



Here, the £ 2 -dependent terms are kept explicitly in order to get the transverse momentum 
dependence of the B-meson WFs. 

By using the equation of motion for the light quark, q lp— 0, and HQET equation of 
motion for the heavy quark, v ■ Dh v = 0, one can obtain four independent differential 
equations, which correlate ^ + (t,z 2 ) and ^-{t, z 2 ) to the 3-particle WFs: 



dV_(t } z 2 ) V+{t,z 2 ) -§_(t,z 2 ) z 2 d 



Of 



d^ + (t,z 2 
dt 



d^ + (t,z 2 
dt 



t dz 2 



$ + (t,z 2 )-*4t,z 2 )} 



udu 



z 2 .~. 



z ~ 



2t[V A (t, u, z 1 ) - V v (t, u, z 2 )} + 3jY A (t, u, z 1 ) - jX A (t, u, z 1 ) (3) 



a*_(t, z 2 ) * + (t, z 2 ) - ^(m 2 ) + 4t _a*+(w, z 2 ) 



dt 



t 



dz 2 



1 

2t 



udu 2t[^ A (t, u, z 2 ) + 2 #v-(t, u, z 2 ) + X A (t, u, z 2 )], 
f z 2 ) - z 2 )! + iAtt+(t, + 2t- a * +( *' ^ 



<9z 2 



rfu (u — 1) 



z 2 ~. 



u, z 2 ) + X A (t, it, z 2 )} - it, z 2 ) + Z —Y A {t, it, z 2 ) 



(4) 



(5) 



and 

* + (M 2 ) *_(M 2 ) 



<9t 



dt 



ik- - 
t 



^+(t,z 2 ) -^_(t,z 2 ) +2t 



dty + (t,z 2 ) 



<9z 2 



<9,2 2 



2i / du(u- 1) $ A (t,u,z 2 ) + Y A (t, u, z )\ 



(6) 



lV '^\^h h \B^^ 1S ^ ne USU& 1 "effective mass" of the B meson in the 



where A = M — 

HQET. By taking the LC limit, one may easily find that our present results agree with 
those in Ref . [12j . Making the Fourier transformation, by virtue of the formulae given in the 
appendix, we obtain, 



UJ 



d^_(cu,z 2 ) _ 2 



du 
du 



t n d 3 

^- A) aJ + 2 



d^ + (LU,z 2 ) d^Jtu,z 2 



dz 2 



dz 2 



+ ^ + (u,z 2 ) = I(uj,z 



^ + {u,z 2 )-^„{u,z 2 ))+A 



d 2 dV+(uj,z 2 



duo 2 dz 2 
d 2 d^ + (co,z 2 ) 



^ + ^z 2 )--^^z 2 ) + 2 duj2 



(7) 

J(uj,z 2 ), (8) 
M(u, z 2 ) + N(u, z 2 )(9) 
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and 



rj2 

A[*+{u>, z 2 ) - z 2 )} + 2— — [*+(«;, z 2 ) + ^ 2 )] = L(w, z 2 ) , (10) 



duodz 2 



where the 3-particle source terms are 



d 



u> roo 



d£ d 



I(u,z 2 ) = 2— dp ^(p,£,^)-*v(P,^ 



-z 2 / dp 



00 d£ ( u — p~~ 



A,-p £ V £ 



[X A (M,^ 2 ) - SY^,^ 2 



J(cu,z 2 ) 
L(c;,2 2 ) 



oi roo 



d£ d 



do; Jo 



-2 f dp f° Ati[Y A ( p ,Z, z z)-2* v (p,t,z 2 )-X A (p,Z,z 2 )] 
Jo juj-d £ at 



d r 00 d£, 



+2— f dpi f [* A (p,£,z 2 )+Y A (p,£,z 2 )] 
duo Jo Jui-p £ 



d 



u roo 



du Jo 
du 2 



d£ d 



M(u,z 2 ) = — dp -± — y A (p,£„z 2 )+X A (p,Z,z 



r dp r ^ A (p,^ z 2 )+x A (p,^z 2 )} 

Jo Jui-p 4 



(11) 
(12) 

(13) 

(14) 



and 



cj roo 



di 



[Y A (p,£,,z 2 )-X A (p,£,z 2 )} 



rui roo 

/ dp / -f[Y A {p,^z 2 )-X A {p,^z 2 )] . 
Jo Ju-p 4 



(15) 



Secondly, by using Eq.(J2J) and the gluon equation of motion, D^G a v {z) = 0, whose source 
term that induces even higher Fock state' contribution is neglected in this work, one can 
obtain two more independent relations among the 3-particle WFs: 

d 



3V A (t, u, z 2 ) + AY A (t, u, z 2 ) - X A (t, u, z 2 



V A (t,u,z 2 ) + 2Y A {t,u,z 2 



1 r.~. 
T 2 



X A (t, u, z 2 )\ + z 2 
+2^$ A (t,u,z 2 ) +Y A {t } u,z 2 )] 



Y A (t,u,z 2 )-X A (t,u,z 2 ) 




ld_ 

tdt 



Y A (t,u,z 2 )-X A (t,u,z 2 ) 



(16) 



and 



j[X A (t,u,z 2 )-Y A (t,u,z 2 )]+2j^[X A (t,u,z 2 )-Y A (t,u,z 2 )} 



-2t 



_d_ 

dz 2 



V A (t, u, z 2 ) + X A (t, u, z 2 )\ - — [V A (t, u, z 2 ) + Y A (t, u, z 2 )\ = . (17) 



Taking the LC limit in Ea. lfTKj) . we get 

m A (t,u)+4Y A (t,u)-X A (t,u) 

with F(t,u) = F(t, u, z 2 )\ z 2^ (F = {^y,^,^^}). By doing the Fourier transformation 
and exploiting the boundary conditions, F(p, £)Ip-^oo,§^oo — ► (F=^ A , X A , Y A ), we obtain 
a relation among the double moments of the 3-particle DAs, 

3 [V A (p, fl]} + 4 [Y^p, 0)1 = [X A (p, 0)1 , (19) 

where [F]^ are double moments of the 3-particle distributions, 

[f% = dp # // j e 1 /••(/'. , (F = {^y, Xl}) . 

J Jo Jo 



+ t 



* A (t,u) + 2Y A (t,u)-X A (t,u) 



(18) 



III. APPROXIMATE SOLUTION FOR THE B-MESON WFS WITH 3-PARTICLE 
FOCK STATES 

In the following, we shall give an approximate solution for the B-meson WFs with 3- 
particle Fock states' contributions by solving the differential equations as shown in Sec. II. 
Before doing this, as a basis and to be self-consistent, we first recollect the results in the 
WW approximation (i.e. I(u) = J{ui) = L(u) = M(u) = N(u) = 0), then, derive several 
constraints for the B-meson DAs <f>±(u)), where the B-meson DAs can be obtained by taking 
the LC limit of the B-meson WFs, i.e. <f>±(w) = lim 2 2^ ^±(^ ; z 2 ). 



A. B-meson WFs ^>^± w {uj, z 2 ) in the WW approximation 

When ignoring the 3-particle Fock states' contributions, i.e. setting I(uj) = J{uj) = 
L(u) = M(u) = N(u) = 0, one can readily obtain the B-meson WW-type WFs ^ w (uj, z 2 ). 
In this case, the two DAs (f)± w (u) take the formP. Il^: 

4>¥ W (") = ^{ U )9{7K - «,); r W ^) = ^^e{u)6(2A - w ), (20) 

and the transverse part, x (u, z 2 ), is a zero-th normal Bessel function. 8(x) is the usual 
unit step function, which equals to for x < and 1 for x > 0. Taking the Fourier 
transformation, ^f±(ou, kjj = / d 2 z± exp(— ik.±-z±)ifj±(u, z 2 )/(2tt) 2 , the normalized B-meson 
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WFs in the momentum space read as 



±ww / \ 

^ w (u, kj.) = ^-^5 fki - lo(2A - u>)) (21) 

7T \ / 



,ww [u] 



and 

k ± ) = ^-^5 (k 2 ± - u(2A - uj) , (22) 

whose k^-dependence correlates to the cu-dependence via a 5-function. 

Eqs. ()21l22|) show that the WFs' dependence on transverse and longitudinal momenta is 
strongly correlated through a combined variable k^/[a;(2A — a;)]. Similar transverse mo- 
mentum behavior for the meson has been discussed in Ref.j^ by transforming the usual 
equal-time wavefunction to its LC form, and has been derived in Ref.j^] by adopting the 
dispersion relations and the quark-hadron duality In these two references, the authors 
stated that the fcT-dependence of the meson's wavefunction depends on the off-shell energy 
of the valence quarks, i.e. ~ k^/x(l — x), where x is the momentum fraction carried by the 
corresponding valence quark. 

B. Some constraints on the B-meson DAs 4>±(lo) 

The B-meson DAs can be obtained by taking the LC limit in equations for the B-meson 
WFs. In the LC limit, Eq.(|7|) is simplified as, 

+0+M=/M (23) 

dui 

with I(u) = I (m, z 2 )\ z 2^ . Similarly, from Eqs. (j4l5J) . we have 

at at Jo 

which leads to 

(u - 2A) </>+(w) + u<j)_(u) = K(uj) (24) 

with 

The solution of the B-meson DAs can be conveniently decomposed into two pieces as 

^)=0rv)+^H , (25) 
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where <f)± w (u) are the DAs in the WW approximation and <j)±\u}) denote what induced by 
the 3-particle source terms I (00) and K(uS). From Eqs. ()23l24|) . the solution for (fy± can be 
obtained straightforwardly, and reads: 



^ ^ = 2A ^ ai u-2A ' ^ ^ 



2A - uj 



(26) 



2A u 

Here, u > 0, a\ and a 2 are integration parameters that satisfy the relation ai + a 2 = 1. 
Note tha » he result in Ref.Q . o* a specific choice of „ = and a2 = 1 in the general 
solution Eq.(j5nj). The function Q{uj) is expressed as follows: 



G(u) = 6{2A - u) U dp 

_ POO 

- 9(u-2A) / dp 



M(p) K(p) 
+ ay 



2A-p 1 (2A-p) 2 

M(p) , K{p) 
ay 



2A J 



2A-p \2A-pf 

with M(p) = I(p) + f i — a2 lp) K{p)- One can easily check that / °° <ko(fy±{uS) = 0, so the 
total DAs are normalized 1 , i.e. Jq°° duxf>±(u)) = 1. 

As usual, we adopt the Mellin moments of 4>±{uj), which take the following form (n = 



0,1,2,---), 



(u n )± = / duj uJ n (j) ± {uo) 
Jo 



duu n (j)^ w {uj)+ I du uj n (j) { l\uj) 
~ Jo 



(28) 



With the help of the formulae given in appendix B, we have 



/, n\WW 
\ W / + 



n\WW 



n + 2 



(2A) n , {u n )1 



(n + l)(n + 2) 



(2A) r 



(29) 



and 



n\(ff) 



(n + l-i)|±_l + lW 



+ (n + 2-i) [X A f~' + (n + 3 - i)^— [^ v ]p 



n + l 



n + 1 , = ! V J / J + 1 



m-1 



(30) 



(31) 



1 Strictly, it is not true 0,0, however since the tail of the LCDAs are proportional to a s , this could be 
accepted as a "tree level" statement. 
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where = — j)!]. The above results for (u> n }± have nothing to do with the free 

parameters a\ and a 2 , due to the fact that a\ + a 2 = 1, and hence they are in agreement 
with the ones obtained in Ref. [12]. It is obvious that the moments of DAs have no relation 
to the 3-particle distribution Y A (p, £)• 

At the present, one knows little about the magnitudes of the B-meson DAs' moments. 
In Ref. js^l , the second moments of the B-meson DAs are estimated by relating them to the 
matrix elements of certain local operators and by calculating these matrix elements from 
the sum rules in HQET, i.e. 

<^ 2 >+ = 2A 2 + 2 -X\ + \x 2 H , <^ 2 >_ = ^A 2 + l -X 2 H . (32) 

Here A^ and Xh parameterize the matrix elements of chromoelectric and chromomagnetic 
fields in the .B-meson rest frame, 

(0\qgE ■ a l5 h v \B(p = 0)> = f B MX% (33) 

and 

(0\qgH ■ cr l5 h v \B(p = 0)> = if B MX 2 H , (34) 
with E l = G°\ H* = -\e ijk G jk , and a = 7°7. From Eqs.Q28j|32J, we obtain 

2 + \ [Xa]\ + = \X% + \xl , [H/ v }\ + \ [X A ]\ = l -Xl . (35) 

Together with Eq. ffEJj) . it leads to 

Wi = |A|, M; = ^ + Ai), [Y a \\ = [X A }{ = -?A|. (36) 



The above results are different from those in Refs. [21, ll2j , where the contributions from [X A ]\ 
and [Y A ]\ have not been taken into consideration 2 . 

As a summary, one may observe that <p±(u) should satisfy the following conditions: 

<fi-(uj)duj =1, / 4> + {uo)duj = 1, (37) 



2 _ r 4 _ 

oj(j)^{uj)duj = -A, / u(j) + (u)du = -A, (38) 

O J Q 

J u 2 ^{u)duj = h 2 + h 2 Hl J uj'M^doo = 2A 2 + 1\% + l -X 2 H) (39) 



2 According to the QCD sum rule analysis [3, as a rough estimation, the contributions to the B-meson WF 
from X A and Y A are at least suppressed by inverse power of the Borel parameter. Here we keep both of 
them for a more complete estimation of the B-meson WFs. 



where according to the QCD sum ru 
Further more, as discussed in Refs. 



e analysis A| = 0.11 ±0.06 and A^ = 0.18 ±0.07. 

DO 

hould satisfy 

A = f — = -L. (X B = 460 ± 160Mey) (40) 

J UJ Aft 

C. Approximate solution for the B-meson WFs ^±(lo,z 2 ) including 3-particle Fock 
states 

To solve the B-meson WFs ty±(u, z 2 ) including 3-particle Fock states, one need to 
know some more details on the properties of the 3-particle WFs \l/y(p, £, z 2 ), ^a(p,(,,z 2 ), 
X A (p, z 2 ) and Y A (p, £,z 2 ), i.e. the transverse momentum dependence of these WFs. In 
the following, we will take three assumptions so as to provide an approximate solution for 
the B-meson WFs ^±{uj,z 2 ) from Eas.pTffl fToTTTI) : 

I) Based on the B-meson WFs in the WW approximation 3,13' we assume that *ff+(u, z 2 ) 
and \fr-(u, z 2 ) have the same transverse momentum dependence, i.e. 

l ± [ W ,z 2 ]=^Mxk^], (41) 

and all the 3-particle WFs also have the same transverse momentum dependence 

X ih) (p^,z 2 ), 

* A (p,t,z 2 ) = * A (p,t)xM(p,ty) ,^v(p^,z 2 ) = MP,0x (h \p,C,z 2 ) (42) 

and 

X A (p,£,z 2 ) =X A (p,0x (h) (p,^z 2 ) ,Y A (p,Z,z 2 )=Y A (p,Ox (h Xp^,z 2 ), (43) 
with the boundary condition lim z 2^ X {P-> £j z 2 ) = 1- 

II) Since the main features of the 3-particle DAs are determined by its first several moments 
(the higher moments will be suppressed by A# or \h accordingly), we assume that the 
relation Eq. (jl9jl among the first non-zero double moments of the 3-particle DAs can be 
extended to be a relation among the 3-particle DAs, i.e. 

Wg WL»M , ( 44) 



10 



71 



III) We adopt a naive model [7( for the difference between \Ev(p, and 9 A (p, £) 3 i 

VvM - *a(p,£) = ^#pe 2 exp . (45) 

With the help of Eqs. (jlfi|17|) and the assumptions (I, II), one can obtain two relations 
among the 3-particle WFs: 

y a ( p , e, z 2 ) = -* A ( P , e, z 2 ) , x A ( P , e, ^ 2 ) = -* a {p, e, A (46) 



1. The transverse momentum dependence of ^±(u>, z 2 ) 

Applying Eq. lfffij) to Eqs.ffFi rroj) . we obtain 
d^Ju.z 2 



UJ- 



(uj -2K)m + (uj lZ 2 ) +u^-(lj,z 2 ) = -L(lu,z 2 ), (48) 
A[^u;,z 2 )-^4u,z 2 ))+2-^[^ + (u;,z 2 )+^4u,z 2 )} = L(^ 2 ) , (49) 

where Eq . (I48j) is obtained from the combination of Eqs. (|8l9|) and the source terms take the 
form, 



'0 -L-p £ \ £ 



^frfJ^VNM,/)] (50) 



and 

L(u,z 2 )=A Tdp r f ^[^(p,£,z 2 )] . (51) 

Substituting Eqs. (j411 142j) into Eq . (|48j) . we obtain a relation between the transverse momen- 
tum distribution of x{ u i z 2 ) an d that of the 3-particle WFs, 



? fo dp c_ p f%[* v (p, qx w (p, e , * 2 )] 
jydpCpf^[*v(p,o] 



- 2 ) = JU r r p ;r - (52) 



It shows that if one knows the 3-particle WF \lV(p, £; -z 2 ) then the exact form of the trans- 
verse momentum distribution x{ u i z 2 ) °f t ne B-meson WFs can be derived; and inversely, a 
constraint on \l/y(p, £, z 2 ) can be obtained as long as one knows the form of xi^i z 2 ). 



3 Even though the first moments of ^>a(p, and ^v(p, m Ref.|]J are different from our's, the difference 
between these two moments are the same for both cases. So we take the same model as the one in Ref.P 
for the difference between 9y(p, £) and \&a(p,£)- 

11 



From Eqs. (j48J) and ()49|) . we have 

_ d 2 



f{oo,z 2 ) = {X-uj)f{u,z 2 ) 



(53) 



dujdz 2 ' 

where f(u,z 2 ) = [^+(uj, z 2 ) + ^-(cj, z 2 )}. Eq.(J53J) shows that the sum of the two WFs 
ty±(u, z 2 ) do not explicitly depend on the 3-particle WFs. Based on the assumption (I), we 
rewrite f(uj, z 2 ) as 



f(u, z 2 ) = + (h(u)z 2 ) = k(u) ■ X (x), 



(54) 



where x = [h(u)z 2 } and the function h(u) is to be determined. Substituting Eq. (j54|) into 
(J53J) . we obtain 

>-A)" 



d 2 x(x) dx(x) 



dx 2 



dx 



1 h(u)K'(u) 
Ii'(u))k(u)) 



2h'(u) 



0. 



(55) 



To ensure the above equation for x{ x ) be tenable to any value of variable u, we set 



1 + 



h!{u)K(u) 





[(w-A)] 


= Cl , 


2/i'(u;) 



C2, 



(56) 



where ci and C2 are two parameters that are independent of variable u. From Eq. (|5l)|) follows 
that 

(57) 



h(u) = ^ '- + c 3 , 



4c 2 

where the value of c 3 is to be determined. 

With the help of Eq. (}56|) . the solution of Eq. (|55|) can be generically expressed as 



X(x)=x(y) = (K 1 r[c 1 }J Cl ^[2y} + K 2 T[2- Cl \J^ Cx [2y]) y 



l-ci 



(cie(0,2)) (58) 



where V is the usual Euler-gamma function, J n is the modified Bessel J-functions, K\ 2 
are undetermined constants and y = ^fo^x (c 2 x > 0). From the boundary condition 
lim 2 2^ x{ u i z 2 ) = 1) we obtain 



lim (KjT^J^y] + ^ 2 r[2 - cj J L _ C1 [2y]) y 1 ^ = 1 
v-*o 



(59) 



Eq . (|59|) leads to K\ = 1 for any value of ci, the value of K 2 is arbitrary for < C\ < 1 and 
equals to for 1 < c\ < 2. Some typical distributions of x{v) are shown in Fig.(JTJ). 

Some discussions about the solution ()58|) for the transverse momentum dependence of 
the B-meson WFs are in the following. 
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y 



FIG. 1: Transverse distributions x(u) with different values of c\ and when < c\ < 1, is fixed 
to be 1. The solid line, the dashed line, the dotted line and the dash-dot line are for c\ = 0.2, 0.8, 
1.2 and 1.8, respectively. 

From the solutions in Eqs. (j571 158|) . one may find that the WFs depend on c<i and C3 only 
in a combined form, so the (C2C3) in practice stands as one free parameter and hereafter, we 
shall always replace (4C2C3) by merely C3 for convenience. 

When C\ = 1 and C3 = 0, we return to the transverse momentum dependence of the 
B-meson wavefunction in the WW approximation. The term including C3 brings some dif- 
ferences to the distribution under the WW approximation, especially the allowed range for 
uj will be broadened for negative value of C3. 

When Ci > 1, one may find the transverse momentum distribution of the B-meson WF 
tends to be a S function as the case in the WW approximation. However, when < C\ < 1, 
the transverse momentum distribution of the B-meson WF will be broadened. This is in 
agreement with the conclusion drawn in Ref.fl^t that the 3-particle contributions might 
considerably broaden the B-meson transverse momentum distribution. To show this point 
more clearly, we take C3 = and transform the transverse part of Eq. (p)Hj) into the momentum 
space. For c\ G (0, 1), we obtain 



1 , r[2-ci]sin[7rci] T A r[2-cj 



^ k±) = - [w^¥ + — — K2 



{{2K-u)u) 



1 



(2A-oj)u> 

where = |kj_|. One may easily find that x(^)kj_) satisfies the normalization condition, 
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FIG. 2: Transverse distributions of x(w, (for = with fixed A = 0.55Gey and uj = 0.5GeV . 
The left diagram is drawn with fixed Ki = 1 and with some different values for ci, i.e. ci = 0.2, 
0.4, 0.6 and 0.8, which are shown in dash-dot, dotted, dashed and solid lines respectively. The 
right diagram is drawn with fixed c\ = 0.6 and different values for K2, i.e. K2 = 0, 5.0, 10.0 and 
15.0, which are shown in solid, dashed, dotted and dash-dot lines respectively. 

/ d 2 kj_x(u, kj_) = 1. The transverse distributions of x(u,k±) with fixed A = 0.55 GeV and 
uj = 0.5 GeV are shown in Fig.Q. The left diagram of Fig.(J2J) is drawn with fixed K% = 1 
and with some different choices for the magnitude of ci, i.e. C\ = 0.2, 0.4, 0.6 and 0.8. The 
right diagram of Fig.(j2J is drawn with fixed c±(= 0.6), and varying K2, i.e. K2 = 0, 5.0, 
10.0 and 15.0. From Fig.(j2J, one may find that with the decreasing of ci, or increasing of 
K 2 , the transverse momentum distributions become broader and broader. 

From Fig.(J2J), one may observe that there is a dip in the transverse momentum distri- 



butions around fcy = y (2A — u)u. It comes from the denominator in Eq. (j60|) . Practically, 

such a dip will not make any problem, due to the fact that we always need the integrated 
results, which will be shown in Sec. IV. By taking the negative value of C3, the allowed 
range of uj will be broadened and it will make a suppression to the singularity in Eq.(|6()jl. 
When summing up all the Fock states' contributions and taking the RG evolution effects 
into consideration, one may expect a further suppression to such singularity in the resultant 
transverse momentum distributions [2j|. 
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2. The distribution functions <f>±(u)) with 3-particle Fock states 



The solutions for the distributions 4>±{oj) with 3-particle Fock states have been given 
in Eqs. (j2(JI25l26|) . However, since the solutions for cp±(uj) (as shown in Eq.(|26|l) involves 
the unknown 3-particle WFs, it can not be used directly. In the following, we shall make 
an attempt to provide more convenient expressions for <f)±(w) under the above mentioned 
assumptions (I, II, III). 

We can derive an expression for the sum of </>±(uj), k(u) = 4>+(u>) + <f>-(u), from 
Eqs. (l56l57J) . which can be expanded in a more convenient form as 

~2A(ci - 1) 



k(u) = c 4 x exp 



-UJ 



x (1 + j3u) 



(61) 



where C4 is an overall normalization factor and we have implicitly taken (— C3) > 0, which is 
reasonable, since the minus sign indicates that the 3-particle WFs will broaden the allowed 
range of u in comparison to that in the WW approximation. Here, f3 is a new phenomeno- 
logical parameter, which stands for the summed effects of other expansion terms 4 . Further 
more, the 3-particle source term can be simplified with the help of Eqs. (|ll|45l46|) as 



\ 2 



— ^(6A 2 - QluA + oj 2 ) exp (-% 
18A 6 V V A 



(62) 



And then the solution for <f)±(u)) can be obtained from Eqs. 



4>+(uj) = C4j3u x exp 



2A(ci 



-u 



-C3 



+ C4,w[/3 + 



2A(ci - 1) 
(-<*) 



Ei 



2A(ci - 1) 



-oj 



-c 3 J 



X 2 



\ 2 



18A 5 



-u exp 



(-u + 5A) + - =- H u Ei — - 
A/ ' 3A 4 V A 



(63) 



and 



)- lu) 



C4 x exp 



2A(ci - 1) 



-oj 



+ 



\ 2 



X 2 



18A 5 



a; exp 



[-uj + 5A) + 



X 2 



3A 4 



x 2 



2A(d - 1) 

UJ 

A 



(64) 



where is an undetermined parameter and the exponential integral function Ei(z) 
— e~ t /tdt. All the terms in the big parenthesis come from the source term I (00). 



4 those terms in higher power series of ui can be summed up, since when u> is big, their contributions are 
suppressed by the overall exponential factor. 
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Some discussions about the solutions ()63I64|) for 4>±{uj) are in the following. 



2A(ci-l) 
(-cs) 



so one can 



Eqs. (j63l64j) show that c\ and C3 are always in a combined form as 
only get the combined results for them. When u — > 0, we have 0+(O) = and 0-(O) = C4. 
Under the condition C3 < 0, one may observe that c\ should be less than 1 to ensure that 
4>±(uj) is normalizable, which shows that the transverse momentum dependence of the B- 
meson wavefunction is broadened due to the introduction of 3-particle wavefunctions. 

From Eqs. (j23l 161(1 . we obtain 

lo^{lo) - 2 4>-{p)dp = f[I{p) - K{ P )]dp. 
Jo Jo 

Numerically, due to the fact that \ 2 E ~ \ 2 Hl Jq I(p)dp « Jq n(p)dp for V oj > 0, and then 
one can safely set = 0, or equivalently X% ~ \ 2 H . 

Under the approximation \ 2 E ~ \ 2 H ~ 2A 2 /3, a solution for <f>±(u) can be directly obtained 
by substituting Eqs. (|63|64jl into constraints ()37p38p39() . i.e. 

0+( w ) = -^exp (-— ) , <j>-{u) = — exp (-— V (65) 



2A(ci-1) 
(-cs) 



_l_ 



where = 2A/3. The undetermined parameters take the following values, 
(3 = Ki = and C4 = — . Such solution for 4>±(uj) also satisfies the constraint (|1D)) , and 
it agrees well with the model for (j)±(uj) raised by Ref.J^, where the same approximation 
\ 2 E ~ \ 2 H ~ 2A 2 /3 is adopted in their QCD sum rule analysis. 



IV. A MODEL FOR THE B-MESON WFS AND ITS PHENOMENOLOGICAL 
CONSEQUENCES 

In the above section, we have derived an approximate expression for the B-meson WFs 
under the assumptions (I, II, III), in which the 3-particle Fock states' contributions are in- 
cluded. The transverse momentum dependence of the B-meson WFs is shown in Eq.([58|) 
and the corresponding DAs are shown in Eqs. ()63l64() . 

For the transverse momentum dependence of the 3-particle WFs, our results indicate 
that when the value of C\ is within the range of (0, 1), it may be expanded to a hyperbola- 
like curve as shown in Figs. (|ll2() . rather than a simple (^-function as is the case of WW 
approximation. Our solution for 4>±{oj) favors c\ < 1 under the condition that C3 < 0, 
which is reasonable since it means that the introduction of 3-particle wavefunctions shall 
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broaden the meson's longitudinal and transverse distributions. This is in agreement with 
the conclusion drawn in Ref.Q], where it has been argued that the 3-particle contributions 
might considerably broaden the B-meson transverse momentum distribution. 

The solutions for the B-meson wavefunction in Eqs. (j58l63l64j) are somewhat complex. 
Based on the discussions in Sec. Ill, we propose a simple model for the B-meson wavefunction 
with 3-particle Fock states' contributions in the following. For convenience, we write the 
two normalized B-meson wavefunctions in the compact parameter 6-space (useful for the 
fey- factorization approach [ill]): 



*+(u,b) = -2 exp (-— J {T[5]Js-i[k} + (1 - S)T[2 - 5]J^ s [k\) ^-J (66) 
and 

*-(w,6) = ^exp (--^-) (T[5]Js-i[^] + (1 - 5)T[2 - 5]J^ s [k]) [^j^ , (67) 



with uj = 2A/3, k = 6(2A — uj)yuj(2A — uj)b and 5 is in the range of (0, 1). In the above 
model, X% ~ ~ 2A 2 /3 is adopted, and to short the uncertainties of the model as much as 
possible, we only take two main phenomenological parameters A and 5 into the definition. 
Here for the transverse momentum dependence part, the range of u is fixed within the 
range of (0, 2A) 5 . When 8 — > 1, the transverse momentum dependence of the B-meson 
wavefunction returns to that of the B-meson wavefunction in the WW approximation. In 
the above definition, the transverse momentum dependence of the B-meson wavefunction 
is still the like-function of the off-shell energy of the valence quarks and keeps the main 
features caused by the 3-particle Fock states, i.e. shall broaden the transverse momentum 
dependence under the WW approximation to a certain degree. 

In order to see the phenomenological influence of the 3-particle Fock states' contributions, 
we recalculate the B — > 7r transition form factor within the /c^-factorization approach and 
show how A and 5 affect the final results. A consistent a naly sis of the B — ► ir transition 
form factor within its physical range has been given in Ref . (Tjj by taking the WW B-meson 
wavefunctions defined in Eqs. (j21l22j) . There we shall adopt the same method as that of 



Ref. 



to do the calculations and to short the paper, we shall only list the results, the 



interested reader may refer to Ref. |l7| for more details of the calculation technology. 



This can be understood by checking the general solution of 4>±(u>) in Eas. (|63lfi4Jl . i.e. the absolute value 
of C3 can not be so big in order for 4>±{lo) to satisfy the constraints (|37I38I39I40|I under the condition of 
ci < 1. 
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FIG. 3: B-meson transition form factors F^(Q 2 ) (Left) and F^(Q 2 ) (Right) with the B-meson 
wavefunction constructed in Eg . (|66l67f) . The upper (lower) edge of the shaded band is for 5 = 0.30 
(<5 = 0.25), and the dash-dot line is for 5 = 0.27. For comparison, the dotted line and the dashed 
line are for the fitted QCD light cone sum rule results is| 
wavefunction [17t |. respectively. A = 0.55GeV. 



and the result derived with WW B-meson 



Naively, the main property of the B — ► 7r transition form factor is determined by the 
first inverse moment of and one may find from Eq.fP that (h^ ] / A^ w ) = (A - 

^b)/^b ~ 0.19, where Aq 9 ^ is the first inverse moment of (u;) and h^ w is that of (f^ w (uo). 
This shows that the 3-particle wavefunctions might be small. In Ref. j^J, by studying the 
B — > ^Iv decay within the perturbative QCD approach, the authors also claims a small 3- 
particle contributions. More explicitly, in Ref . [2J| , the 3-particle contributions are estimated 
by attaching an extra gluon to the internal off-shell quark line, and then (l/m^) power 
suppression is readily induced. In the following, we shall study the uncertainties caused by 
two parameters S and A under the condition that the 3-particle contribution is less than 
±20% of that of the WW case, and at the same time give the possible range for S and A. 

By taking B-meson wavefunctions as Eqs. (j66pri7|) . we first study the uncertainties of 
B — > 7i transition form factor caused by S with fixed A = 0.55Ge\^ (the center value of A 
determined in Ref. \v\ ) . We show the B — > ix transition form factors F^q(Q 2 ) in Fig.Q. 
Our results show that if the contribution from the 3-particle wavefunction is limited to be 
within ±20% of that of WW wavefunction with Q 2 G (0, ~ lOGeV 2 ), then the value of S 
should be within the region of (0.25, 0.30). 

Next, we study the uncertainties of B — ► tt transition form factor caused by A with fixed 
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FIG. 4: B-meson transition form factors F^(Q 2 ) (Left) and F^(Q 2 ) (Right) with the B-meson 
wavefunction constructed in Ea. l|66l67|) . The upper (lower) edge of the shaded band is for A = 0.52 
(A = 0.58), the dash-dot line is for A = 0.55. For comparison, the dotted line and the dashed line 



18j | and the result derived with WW B-meson 



are for the fitted QCD light cone sum rule results 
wavefunction [Tjj], respectively. 5 = 0.27GeV. 

5 = 0.27 and the results are shown in Fig.Q. Our results show that if the contribution from 
the 3-particle wavefunction is limited to be within ±20% of that of WW wavefunction with 
Q 2 G (0, ~ lOGeV 2 ), then A should be within the region of (0.52, 0.58)GeV. Similarly, one 
may find that if taking a larger value for 5 (e.g. 5 = 0.30), then the range of A should be 
shifted to a bigger interval (e.g. (0.55, 0.61)GeV). 

Figs. (011) show that by taking into account the 3-particle wavef unctions' contributions, 
the B — > 7i transition form factors F+q(Q 2 ) raise slower with the increment of Q 2 than the 
case of WW B-meson wavefunction. And if the contribution from the 3-particle wavefunction 
is limited to be within ±20% of that of WW wavefunction with Q 2 G (0, ~ WGeV 2 ), then 
the possible range of 5 and A are, 5 ~ (0.25, 0.30) and A ~ (0.50GW, 0.60GeV). 

V. SUMMARY 



It had been proved that the B-meson WF is renormalizable after taking into account the 
RG evolution effects Q. a nd the —le featute Q of toe B-nteson DA can be retnoved 
under evolution. Therefore, to keep the kx dependence in both the hard scattering ampli- 
tude and the wavefunction is necessary. It was found that the transverse and longitudinal 
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momentum dependence in the B-meson WF under the WW approximation is correlated 
through a 5-function, S(h 2 ± — uj(2A — u)). In the paper, we show that the transverse mo- 
mentum distribution of the B-meson WF can be broadened to be a hyperbola-like curve by 
including 3-particle Fock state, rather than a simple ^-function. 

The solutions in this paper provide a practical framework for constructing the B-meson 
LC WFs ty±(u), z 2 ) and hence are meaningful for phenomenological applications. And we 
have constructed a new model for the B-meson wavefunction in the compact parameter b- 
space as shown in Eqs. (j66l67J) based on these solutions. There are uncertainties caused by 
two unknown parameters A and S. However, since the B-meson WFs are universal, we can 
determine them by global fitting of the experimental data. By taking B — > 7r transition 
form factor as an example, we show that if the 3-particle wavef unctions' contributions are 
less than 20% of that of the WW case, then one may observe that the preferable values for 
these two parameters are S ~ 0.27 and A ~ 0.55GeU. 

The reasonable inclusion of the 3-particle Fock states in B-meson WFs provides us with 
the chance to make a more precise evaluation on the B meson decays. Further studies on 
the B-meson WFs with higher Fock states and its phenomenological implications are still 
necessary. 
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APPENDIX A: BASIC FORMULAE FOR THE FOURIER TRANSFORMATION 

First, we define F(t,u) = F(t,u, z 2 )\ z 2^ and the Fourier transformation: 



Here uv + and £v + (v + is the '+ '-component of v defined in Eq.Q) denote the LC projection 
of the momentum carried by the light antiquark and the gluon, respectively, and F(u, £) 
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vanishes unless u> > and £ > 0. 
Some useful formulae: 

^L e ^(t,z 2 ) = ^,z 2 ), (A3) 

dt iuJt d^{t,z 2 ) . T . 2 , ,. A . 

— \i± = A4 

-tV<^(t, = (-i)»_*(u;, A (n > 1) (A5) 

^^-^T 1 = (« > 1) (A6) 

where \I> = (^ r + ,^ r _) and \& = (\E f + ,\I / _), respectively. And for the 3-particle distributions, 
we have 

*%*e**9(t,u) = fdpT %*(p,t), (A7) 

Z7T JO Ju)-p C, 

J ^tV**(t, u) = (-iT-^ £ dp jT f 0, (n > 1) (A8) 

e^{t,u) = dp -1 _±U( p ,£), (A9) 



27T JO Joj-p £ \ £ 



27T JO Ju-p £ <9£ 

where \I/ = (^V, ^a, Xa, Ya) and ^ = (&v, ^A, Xa, Ya), respectively. Note here, we have 
implicitly using the following equation, 

roo poo rl ruj roc AC 

\ dp dU du*{p,0-6(u-p-Zu)= dp (A12) 

JO JO JO JO Ju-p £ 

APPENDIX B: MELLIN MOMENTS OF THE DISTRIBUTION AMPLITUDE 



In order to calculate the Mellin moments defined in Eq. (j28|) . it is more convenient to use 
the derivative of Q{uj), i.e. 



d 2A 



dui uj{2K — uj) 

And from Eq. ()23|) . we have the following equation 



a x K{u) a 2 K(u) d 

A a; ) + 7TT 1 a 2 —K{uj) 

2A-oo uj duo y ' 



(Bl) 



(u,»)_ = ^!Z± L_ / u"l( u )du. (B2) 



n + 1 n + 1 jo 
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Substituting En. fjTTJl into Eq.(jB2j), we get 

d 



poo roc (f poo poo pi r) 

/ co n I(uo)duo = 2 / u; n — - / dp/ df / dit— [* A - * v ] • <J(w - p - fit 
JO jo du Jo Jo JO Oq 



du 



yn-l 9 



poo poo pi 

-2n / dp di du( P + tuy-^iVA - ^ 

Jo Jo Jo af 



oo roo 



/o JO 
n-1 



E7 / n — 1 
, ,7 + 1 



n-1 



where the double moments of the 3-particle distributions are defined as 

coo roo 

dp / dfp* -X J 'I'd'. f ) (F = {*v, *A, X A }) 
o Jo 



(B3) 



(B4) 



The moments of 4>+\u) can be written as 



l r°° u n+1 



n + 2 Jo 2A-u 

roo 



+ -Zt + a 2 —K{u) 

2A — u u duo 



du 



+ai 



o 



K{u)du . 



u-2A 

More definitely, with the help of the Eqs. (jllJ12)) . we have 



(B5) 



n\{g) 



u n+1 d 



n + 2 J 2A - udu 



+ 



U 



n+1 f d 



n + 2 7 (2A-u;) 2 [du; 

+2 

r d 



d 

dpd£du—[ty A - * v ] -6(u-p- fit) 



dpdfdwf^A + -^a] • ^(^ — P — f M ) 



du 



2A — u duo 



+2 



5 

dpdfdit— |Wy] • <5(u; — p — fit) 



dpdfdw[\l/yi + X^] • 5 (a? — p — fit) 



9 

dpd^du— [\£v] • 5 (a; — p — fit) 



da; 



du; 



where we have implicitly applied the relation ai + ct2 = 1. To get the final results, the 
following formulae are useful (n > 2), 



do; 



00 roo rl Q 



" 



n-1 



j (n-1 



dp df du-[F].8{u-p-^u) du = nE-^TT ■ W 
./o Jo df fel J + 1 V J / 



m-i 



°° u n — 
duo 



OO /"OO /"l 

<^P / ^f / dit[F] • 5 (a; — p — fit) 
jo jo 

00 roo rl Q 



roo roo roo rl H q T)\ 



do; 
du; 



El"|m 



i=i 



(■00 



du; 2 



/■OO /"OO cl "1 ™ 1 

/ dp df / dit[F] ■ 5(u; — p — fit) du; = n 
jo jo jo J 



n — 1 



in; 



n-1 
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where F = (9 V , V A , X A , Y A ). 

With the help of the above formulae, the value of (u n ) ( f > can be directly derived, as is 
shown in Eq. (J3U)) . One subtle point is that, before doing the integration over u, it is more 
useful to expand 2 ^ 1 _^ , i.e. 

-i( 1+ m + m\(*A)\.., (B6) 



2A — uj uj \ \ uj J \ uj ) \ uj 

where the power of (j^j will be stopped at a particular value, for one may observe that the 
terms with even higher powers contribute zero exactly. And to do the integration in a form 
like 

, /^W (B7) 

(w-2A) 2 duj K ' 

where TC(oj) is a function of u, we can transform it to a more familiar one 

f u n dH(u) J f uu n ~ l dH{u) J f uo n d 2 H{uu) J 

/ i r^7 — T^duo = n / -rr- — r^duj + / — — -r^du . (B8) 

J (uj — 2A) du J (uj — 2A) du J (u - 2A du 2 v ; 
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